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Abstract 

We discuss a Lattice QCD mixed action investigation employing Wil- 
son maximally twisted mass sea and overlap valence fermions. Using 
four values of the lattice spacing, we demonstrate that the overlap 
Dirac operator assumes a point-like locality in the continuum limit. 
We also show that by adopting suitable matching conditions for the 
sea and valence theories a consistent continuum limit for the pion de- 
cay constant and light baryon masses can be obtained. Finally, we 
confront results for sea-valence mixed meson masses and the valence 
scalar correlator with corresponding expressions of chiral perturbation 
theory. This allows us to extract low energy constants of mixed ac- 
tion chiral perturbation which characterize the strength of unitarity 
violations in our mixed action setup. 

PACS numbers: 11.15.Ha, 12.38. Gc 
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1 Introduction 



The discovery that despite the no-go theorem of Ref. [T], an exact chiral sym- 
metry can be established for lattice regulated field theories without violating 
essential field theoretical conditions [2], is certainly one of the major concep- 
tual breakthroughs within lattice field theory in the last years. The funda- 
mental theoretical development for this achievement is the Ginsparg- Wilson 
relation [3]. The theoretical and conceptual advantages of a lattice Dirac 
operator obeying this relation have been pointed out in Refs. [U [5j [6j [7] . A 
practical solution - suitable for numerical simulations - of an operator satisfy- 
ing the Ginsparg- Wilson relation, termed overlap Neuberger- Dirac operator, 
was found in Refs. |SJ [H]. 

Unfortunately, soon after this exciting development in lattice field the- 
ory, it turned out that for numerical simulations overlap fermions are very 
expensive. In fact, the numerical cost is so demanding that many groups 
working in lattice QCD are nowadays still using lattice discretizations based 
on computationally much cheaper Wilson or staggered like fermiong^. Al- 
though with such much simpler lattice fermions, presently, simulations very 
close to or even at the physical value of the pion mass can be performed, see 
e.g. Ref. j 13 j . it would be nevertheless highly desirable to take advantage of 
the conceptually much cleaner overlap fermions, or other Ginsparg- Wilson 
fermions that can be applied in practice, to investigate problems where chiral 
symmetry or topological aspects play a significant role. Given the high nu- 
merical demand of overlap fermions, such simulations are, however, restricted 
so far to the future, when supercomputers of much higher performance than 
available today, in combination with, hopefully, much improved algorithms 
are developed. 

An intermediate step towards the goal to employ overlap fermions is to 
use them only in the valence quark sector, while generating gauge field config- 
urations with much cheaper sea quarks, e.g. Wilson or staggered like quarks. 
Such a situation is called mixed action setup. The basic idea in such a setup 
is to relate the sea and valence quark actions by a suitable matching condi- 
tion, e.g. by demanding that at a given value of the lattice spacing the pion 
mass is the same in both discretizations. With such a condition, in principle, 
the continuum limit of a mixed action setup is well defined and should give 
correct answers for all physical observables of interest. 

However, at any non-zero value of the lattice spacing care has to be taken 
when using mixed actions. Since the sea and valence actions are different, 

1 See, however, the work in a chiral invariant Higgs- Yukawa model, where extensive use 
of the overlap operator is made |10l HT1 [TS] . 
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there are unitarity violations the size of which is unknown a priori and should 
be determined. Also, the eigenvalue spectra of the sea and valence lattice 
Dirac operators are not matched, which can lead to problems especially for 
the zero mode contributions of the chiral invariant overlap Dirac operator. In 
addition, the locality property of the overlap Dirac operator [14j needs to be 
determined in order to see whether hadronic quantities could be affected by 
finite values of the exponential decay rate. All these aspects demand therefore 
a principal investigation of a mixed action before large scale simulations can 
be started. 



In this paper, we want to report on one such basic mixed action analysisj, 
namely using maximally twisted mass fermions [25] in the sea sector and 
overlap fermions in the valence sector. In a previous work [26, 27J, a first 
study of this mixed action setup has been performed and particular difficulties 
with such an approach have been reported, which arise from the fact that 
the exact zero modes of the overlap valence Dirac operator are not matched 
exactly by the sea quark twisted mass Dirac operator. This phenomenon 
appears at typical values of the lattice spacing currently used in simulations, 
i.e. 0.04fm < a < O.lfm. In such a situation, special care must be taken, in 
particular in the way the sea and valence quark theories are matched. 

Here we want to extend the analysis of Ref. j2jj in several directions. The 
first is that we add a new and finer value of the lattice spacing than available 
in Ref. [27J to the analysis of the continuum limit scaling test of the pion 
decay constant. As we will discuss below, the results of the analysis with 
this new value of the lattice spacing confirm the findings of Ref. [27] and 
strengthen the conclusion that there are critical values for the pion mass and 
the physical volume below which simulations can strongly suffer from the 
effects of the mismatch of the zero modes between sea and valence sectors. 
We will also demonstrate that the improved matching condition, suggested 
in Ref. [27J, indeed works well for the continuum limit of the pion decay 
constant and the nucleon and A masses which are new observables added to 
the study of the properties of the mixed action we are considering here. 

As a second aspect, we want to show that the overlap Dirac operator, 
which is only exponentially localized at non-zero lattice spacing [14], reaches 
the expected point-like locality property in the continuum limit. We will 
use four values of the lattice spacing to perform such a continuum limit 
extrapolation and, in particular, compare the localization range of the overlap 

2 Mixed actions have been used by various groups. Here we mention studies using 
domain wall sea and overlap valence fermions |15[ 1TB] . Wilson clover sea and overlap 
valence quarks |17l ITS] . domain wall valence fermions on a staggered sea |19l 1201 I2T] . 
Moreover, another class of mixed actions involves using variants of the sea action in the 
valence sector, such as Osterwalder-Seiler quarks on a twisted mass sea [2211251123] , 
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Dirac operator with hadronic scales, i.e. the pion and the nucleon masses. 
To our knowledge, this is the first unquenched continuum limit study of the 
locality properties of the overlap Dirac operator. 

As a third and also new investigation, we study the meson masses con- 
structed in the mixed action theory and confront them with expressions of 
mixed action chiral perturbation theory. This allows us to compute a num- 
ber of low energy constants of the mixed action effective chiral Lagrangian. 
In addition, we will look at the non-singlet scalar correlator to test for pos- 
sible unitarity violation which are inherent in mixed action simulations, as 
mentioned above. 

The main goal of the paper is to provide a basic and principal investigation 
of a mixed action employing overlap fermions in the valence sector. The 
results of this investigation can and will help any further calculation aiming 
at computing physical quantities in that it provides important information 
on the properties of mixed actions which can help to perform safe simulations 
in the future. 

The outline of the paper is as follows. Sec. 2 discusses some theoretical 
aspects of our setup and gives parameters of our lattices. In Sec. 3, we look at 
the continuum limit behaviour of the locality of the overlap Dirac operator. 
Sec. 4 reports the results of our new and finest lattice spacing by including 
them in a continuum limit scaling test of the pion decay constant. In Sec. 5, 
we discuss mixed meson masses and unitarity violations in the non-singlet 
scalar correlator. In Sec. 6, we employ our mixed action setup in the baryon 
sector and perform a continuum limit scaling test for nucleon and A masses. 
Sec. 7 concludes and summarizes our main findings. 

2 Simulation setup 
2.1 Overlap Fermions 

The fundamental Nielsen-Ninomiya theorem [UEE], formulated in 1981, ap- 
parently excluded the possibility of simulating chiral fermions on the lattice, 
without violating essential properties, such as the absence of doubler modes 
or the locality of the theory. However, already in 1982 it was shown by 
Ginsparg and Wilson j3] that a remnant of chiral symmetry is present on the 
lattice without the doublers, if the corresponding lattice Dirac operator D 
obeys an equation now called the Ginsparg- Wilson relation: 



For many years, though, no practical solution to this equation had been 
known. The situation changed in the second part of the 1990s when the 



75j D + D75 = aD lb D. 
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Ginsparg- Wilson relation was reconsidered and its conceptional advantages 
pointed out [U El El IS]- Moreover, in Ref. [2] the key observation was made 
that any lattice Dirac operator satisfying the Ginsparg- Wilson relation leads 
to an exact chiral symmetry at non-vanishing value of the lattice spacing. 
Neuberger found a particular closed form of a lattice Dirac operator D - 
called overlap Dirac operator - that obeys the Ginsparg- Wilson relation and 
that can be employed in practical simulations. The massless overlap Dirac 
operator is given byH [H |9]: 

D ov (0) = i(l-A(AU)- 1 / 2 ). (2) 

In the kernel operator A we choose the standard Wilson-Dirac operator: 

A= l + s-aD WilBon (0), (3) 

where s is a parameter which satisfies \s\ < 1 and can be tuned to optimize 
locality properties [14| (see Sec. [3] for a detailed test of locality in our setup). 
The Wilson-Dirac operator is defined by: 

Aviison(™ ) = - (ta»(V; + V„) - aV^V^) + m , (4) 

where mo is the bare Wilson quark mass and V M (V* ) are the forward (back- 
ward) covariant derivatives. The massive overlap Dirac operator with the 
bare overlap quark mass m ov reads: 

A>v(m ov ) = (l - ^y 1 ) Aw(0) + m ov . (5) 

The most important property of overlap fermions is the fact that chiral 
symmetry in a lattice modified form can be established with very important 
consequences, among them the absence of 0(a) lattice artefacts. 

From the technical side, an important feature of overlap fermions is that 
the explicit construction of the overlap Dirac operator involves an approxima- 
tion (in our case the Chebyshev polynomial approximation) of the operator 
(A^ A)~ 1 / 2 that needs to be realized up to machine precision, in order to 
consider the lattice chiral symmetry to be exact. This leads to a large com- 
putational cost of overlap fermions, around two orders of magnitude larger 
than of e.g. Wilson twisted mass fermions. For a comparison of overlap and 
Wilson twisted mass fermions in terms of computing cost, we refer to |29j . 

Further problems occur in simulations with dynamical overlap fermions. 
The overlap Dirac operator develops discontinuities when changing topologi- 
cal sectors and the proposed solutions to alleviate this problem [3T> 1 13"T ] I3"2" | I3"5] 

3 For an early review of the overlap fermions formalism, see e.g. Ref. [7]. 
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lead to a further growth of the computational cost. The discontinuity prob- 
lem can also be avoided by fixing topology [34, 35J, which, however, still 
needs very computer time expensive simulations and introduces additional 
0(1/ V) finite volume effects. The mixed action approach, which uses compu- 
tationally cheap sea quark simulations on which then valence overlap quarks 
are evaluated, is another way of avoiding the large computational costs of 
dynamical overlap simulations. This is the setup considered in this paper. Of 
course, in this approach exact chiral symmetry is only realized in the valence 
sector. 

Since the pseudoscalar decay constant will play an important role in the 
following, we discuss this observable here. For overlap fermions, the pseu- 
doscalar decay constant f° v can be extracted using the PCAC relation and 
does not require the computation of any renormalization constants (see e.g. 
Ref. |3BJ): 

/r = ^l(0|P|vr) ov |, (6) 

where m™ is the mass of the charged pseudoscalar meson and |(0|P|7r) ov | the 
matrix element of the pseudoscalar current, both extracted from the two- 
point pseudoscalar correlation function Cpp(t), built of two mass-degenerate 
overlap valence quarks. 



2.2 Wilson Twisted Mass Fermions 

Wilson Twisted Mass (tm) fermions [37] were originally introduced to deal 
with the problem of unphysically small eigenvalues (zero modes) of the Wilson- 
Dirac operator. As pointed out in Ref. [25] an essential advantage of this 
formulation of lattice QCD is also the possibility to obtain automatic 0(a)- 
improvement, by tuning just one parameter - the bare Wilson quark mass 
to its critical value. This property has been confirmed in detailed contin- 
uum limit scaling studies in the quenched approximation [38] and with two 
dynamical quarks [39]. What is more, the twisted mass discretization can 
reduce the effects of explicit chiral symmetry breaking by the Wilson term 
in the renormalization process. In fact, the problem of operators belonging 
to different chiral representations can be suppressed or even avoided in this 
formulation. Among the disadvantages of this formulation are the explicit 
breaking of parity and isospin symmetry, being, however, 0(a 2 ) cut-off ef- 
fects [251 EH El 1121 [39] , which were observed to be substantial only in the 
neutral pion mass. 

Twisted mass fermions are obtained by adding a chirally rotated mass 
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term to the Wilson-Dirac operator in Eq. (j3J) in the following way: 

Am = Avilson(rao) + Witt , (7) 

where fi is the twisted mass parameter and r 3 is the third Pauli matrix acting 
in flavour space. 

A comprehensive simulation programme has been undertaken by the Eu- 
ropean Twisted Mass (ETM) collaboration, including simulations with dy- 
namical up and down quarks (JVf = 2) [l3j HU [39] and with dynamical up, 
down, strange and charm quarks (iVf = 2 + 1 + 1) |4"5] |4"S] . 

The charged pseudoscalar meson decay constant can be obtained in a 
similar way as in the case of overlap fermions, without the need to compute 
any renormalization factors: 

^ = ^l<0|P|^tJ, (8) 

where m^ m is the pseudoscalar meson mass and | (0| J => |7r =l= ) tm | the matrix el- 
ement of the pseudoscalar current, both extracted from the two-point pseu- 
doscalar correlation function Cpp(t), built of two mass-degenerate tm valence 
quarks. 

2.3 Mixed action setup 

Our mixed action setup consists of two mass-degenerate flavours of Wilson 
twisted mass quarks at maximal twist in the sea sector and overlap valence 
fermions. The parameters of our dynamical ensembles are provided in Tab.[TJ 
The ensembles were generated by the ETM Collaboration [39], using the 
tree-level Symanzik improved gauge action and Nf — 2 flavours of twisted 
mass fermions, tuned to maximal twist by setting the hopping parameter 
k = 1/(8 + 2amo) to its critical value, at which the PCAC quark mass 
vanishes [43j. 

The light-quark ensembles, labeled in Tab. [1] with a subscript £, cor- 
respond to a fixed physical situation with roughly fixed pseudoscalar me- 
son mass m^ro ~ 0.8 (where r is the Sommer parameter |47j ) and lat- 
tice size L w 1.3 fm. The non-perturbatively renormalized |48j quark mass 
/i^ s (/i = 2 GeV) ~ 20MeV, which in infinite volume gives a pseudoscalar 
meson mass of around 300 MeV. The chirally extrapolated values of the 
Sommer parameter for our ensembles are r^/a = 5.25(2) at = 3.90, 
r / a = 6.61(2) at p = 4.05 and r /a = 8.33(5) at /? = 4.20, which cor- 
responds to lattice spacings a rs 0.079 fm, a ~ 0.063 fm and a 0.051 fm, 
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respectively [31]. At (3 = 4.35, we only have one sea quark mass available, 
which yields r /a = 9.82(4) and a lattice spacing a ~ 0.042 fm. 

For our investigations of locality (Sec. |3j, we use ensembles with L/a = 32 
at: = 3.9, afi = 0.004 (ensemble B i>32 ), (3 = 4.05, afi = 0.003 (C i>32 ), 
(3 = 4.2, afi = 0.0065 (with a heavier quark mass, ensemble -D/^32) and 
(3 = 4.35, afi = 0.00175 (Eg). In addition, we use two additional volumes at 
(3 = 3.9, afi = 0.004 (with linear lattice extent of 0.6 fm (Bg$) and 1.6 fm 
[Bg, 20)), to check finite size effects in the context of locality. 

We have also used one small-volume ensemble with a heavier quark mass, 
labeled as B^. It corresponds to the same lattice size and lattice spacing as 
Bg, but the renormalized quark mass is fi^ s (fi = 2 GeV) ~ 40MeV, which 
corresponds to m^r^ w 1.0 and a pseudoscalar meson mass m n ~ 450 MeV 
in infinite volume. 

Finally, to explicitly check whether one is safe against the effects of zero 
modes of the overlap operator, discussed in Ref. |27| . we have also considered 
an ensemble B s that corresponds to a physical extent of the lattice of L rs 2 fm 
and a rather heavy renormalized quark mass of fi^ s (fi = 2 GeV) ~ 45 MeV, 
which yields m n r « 1.0 and a pseudoscalar meson mass of « 480 MeV 
in infinite volume. For this ensemble, m n L m 4.7 and hence the effects of 
zero modes should be very much suppressed |27| . 

Our setup for the overlap fermion valence sector has been described in 
detail in Ref. |49] . Here we only briefly mention that we performed one 
iteration of HYP smearing [50J of gauge links to reduce the condition number 
of A* A. To assure the best locality properties, we set s = in Eq. ([3]) (for 
details, see the next section). 

We have produced all-to-all overlap and twisted mass propagators for a 
wide range of quark masses, which allowed us to compute three types of 
correlation functions: sea-sea (unitary tm), valence- valence (overlap) and 
valence-sea (mixed overlap-tm). 

3 Locality of the overlap operator 

The aim of this section is a comprehensive investigation of the locality of the 
overlap lattice Dirac operator. We will use ensembles at four values of the 
lattice spacing, which will allow us to perform a continuum limit study of 
the effective decay rate of the overlap operator to see whether indeed point- 
like locality is recovered in the continuum limit. We will also investigate 
the influence of the parameter s in Eq. ([3]), as well as finite size effects, for 
which we use four different physical volumes. Finally, we will perform a 
comparison of the decay rate with hadronic scales of the theory, in particular 
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Label 


(3 


(L/af x T/a 




^crit 


m^ m r 


L/r 


m tra L 


Be 


3.90 


16 3 


x 32 


0.0040 


0.160856 


0.84 


3.0 


2.5 


c t 


4.05 


20 3 


x 40 


0.0030 


0.157010 


0.83 


3.0 


2.4 


Dt 


4.20 


24 3 


x 48 


0.0020 


0.154073 


0.82 


2.9 


2.4 


Ei 


4.35 


32 3 


x 64 


0.00175 


0.151740 


0.74 


3.3 


2.4 


B h 


3.90 


16 3 


x 32 


0.0074 


0.160856 


1.03 


3.0 


3.1 


B s 


3.90 


24 3 


x 48 


0.0085 


0.160856 


1.02 


4.6 


4.7 


Be,s 


3.90 


8 3 


x 16 


0.0040 


0.160856 




1.5 




Be,2o 


3.90 


20 3 


x 40 


0.0040 


0.160856 


0.73 


3.8 


2.8 


Be,32 


3.90 


32 3 


x 64 


0.0040 


0.160856 


0.70 


6.1 


4.3 




4.05 


32 3 


x 64 


0.0030 


0.157010 


0.69 


4.8 


3.3 


B>h,32 


4.20 


32 3 


x 64 


0.0065 


0.154073 


1.10 


3.8 


4.2 



Table 1: The parameters of dynamical Nf — 2 maximally twisted mass en- 
sembles used in this work. We give the label, the values of the inverse cou- 
pling (3, the lattice volume L 3 x T, the twisted mass parameter a/i, the 
critical hopping parameter K cr ; t = 1/(8 + 2am cr i t ), the approximate values of 
the pseudoscalar meson mass and of the lattice size in units of the Sommer 
scale r and the product m^ m L. As demonstrated in Ref. [27], the values of 
m^ m L provide a very good estimator for the the role of the zero modes of the 
overlap operator: they can be large (if m^ m L < 3), potentially large in some 
observables (if 3 < m^ m L < 4) or almost negligibly small (if m^ m L > 4). 

the pseudoscalar and the nucleon masses, in order to achieve a quantitative 
measure for the decay rate relative to physical scales. Such a comparison 
is clearly of importance: as long as the decay rate is smaller than, say, a 
hadron mass, computations with the overlap operator are not very useful. 
If, as expected, the overlap operator converges in the continuum limit to the 
Dirac operator with point-like locality, there must exist a value of the lattice 
spacing, where the decay rate is much larger than the hadronic scale. It is 
therefore one goal of this paper to determine the ratio of the hadronic mass 
and the decay rate and to determine this ratio at the values of the lattice 
spacings used here. 

In this section, we will closely follow the investigation of the locality prop- 
erty of the overlap operator of Ref. [H]. While for small enough values of 
the lattice spacing, the locality of the overlap operator can be analytically 
proven, for coarser values the locality property needs to be determined nu- 
merically. In fact, in Ref. [2] it has been demonstrated that the locality 
property can be maintained to rather coarse values of the lattice spacing of 
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about a w 0.1 fm, which has been pushed to even coarser values of a ~ 0.2 
fm in [51]. Since the values of the lattice spacings employed in this work are 
all below 0.1 fm, we expect that the valence overlap operator is exponentially 
localized for all ensembles considered, an expectation which will be confirmed 
below. 



3.1 Effective decay rate 

It is clear that the overlap Dirac operator (denoted D in this section) is not 
local in a geometrical sense, which would mean that there exists a finite £ 
such that the following requirement is fulfilled: 

Dtf (n) = D ( n , m)*(m), (9) 

m 

D(n, m) = for all points n, m such that ||n — m|| > £, 

for some norm ||-|| to be defined below. However, from the point of view 
of continuum limit behaviour of physical quantities, such geometrical local- 
ity is not necessary. It is sufficient if the Dirac operator D(n,m) decays 
exponentially fast, i.e.: 

||£>(n,m)|| < Ce^ l|n - w|1 (10) 

with some decay rate p (in lattice units, with corresponding physical decay 
rate p = p/a) and normalization factor C. 

For the norm in Eq. (191) , we use the taxi-driver distance for periodic 
boundary conditions: 

IM1 X = y^min{ln M | , - n M |}, < < jV M , (11) 

where iV M is the number of lattice sites in direction p. We also define the 
norm of the operator as the row sum norm: 



4 3 



\\D(n,m)\\ = max . > > \D(n, m)g (12) 

l<a<3 y=1 6=1 

where /i, z/ denote the Lorentz indices and a, 6 the colour indices. For the 
same value of the taxi-driver distance, the operator norm can take several 
values. Hence, we also define the maximum norm: 

IPIL*(d)= m ^x ||D(n-m,0)||. (13) 

\\n—rn\\ 1 =d 

From now on, we will only consider the maximum norm, since it corresponds 
to the most restrictive case, i.e. proving locality using the maximum norm 
implies locality for any other choice of the norm. 
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3.2 Continuum limit analysis 

In this subsection, we will perform a continuum limit scaling test of the decay 
rate of the overlap operator using ensembles B^, CV,32, -D/1,32, Ei of Tab. [TJ 
These ensembles correspond to four lattice spacings, at different physical 
volumes, but keeping constant the lattice size L/a = 32. 

The Dirac operator is local if it fulfills Eq. (ITUj) . We define the effective 
decay rate p e g as follows: 

^ = ^( nm l 7 in +yX ( 14 ) 

Larger values of the decay rate imply better locality properties. 

The full decay property of the overlap operator, including small distances, 
is determined by a sum of exponential decays with corresponding decay rates: 

PLax = £^ e ~*' M \\D\\ max = Ce-'HI . (15) 

i 

Therefore, the desired decay rate is only attained at asymptotically large 
distances, where the effective decay rate p e g will be constant as a function of 
the taxi-driver distance. 

In Fig. [TJ we show the effective decay rate as a function of the taxi-driver 
distance for four different values of the lattice spacing. Our typical statistics 
was around 50 configurations well separated in Monte Carlo time. In order to 
observe the expected plateau-like behaviour of p e g, it is necessary to have a 
large value for L/a. If L/a is too small, the signal is dominated by hyper cubic 
artefacts and the extracted values of p are not reliable. We have studied this 
effect in free- field theory, for different values of L/a and we have concluded 
that L/a = 32 is necessary to observe a stable plateau of the effective decay 
rate (see Sec. 13.51 for details). 

Moreover, from Fig. [T] it is clear that the determination of the effective 
decay rate is not straightforward, i.e. the observed large fluctuations have to 
be taken into account in the analysis as a systematic error. To this end, we 
follow the strategy of Refs. |39 [ 152" ] . As a first step, we fit the effective decay 
rates obtained at different values of the taxi-driver distance to a constant, 
including at least 4 consecutive points. For this, we fit the value of the decay 
rate and the corresponding statistical error using a jackknife procedure. For 
each fit, we also compute the corresponding value of \ 2 and the confidence 
level CL, defined as: 

CL(g = X \ n = d.o.f.) = 1 - 7 (g/2, n/2), (16) 
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Figure 1: Effective decay rate p [fm -1 ] as a function of the taxi driver distance 
[fm] for four different values of the lattice spacing. The blue line is the final fit 
obtained following the method explained in the text. The fitting intervals [fm] 
are indicated in Tab. [3j The error bars are often smaller than the symbol size; 
therefore the large fluctuations induce significant systematic uncertainties in 
the determination of the effective decay rate. 



where d.o.f. denotes the number of degrees of freedom for a given fit and 
7(g/2,rz/2) is the lower incomplete Gamma function. 

In this way, we obtain a large sample of fitted values of the effective decay 
rate from which we can construct weighted distributions using two kinds of 
weights, the confidence level (CL) or a function of \ 2 which we choose as 
exp(x 2 )/d.o.f. - thus damping the influence of bad fits. From this weighted 
distribution, we extract the mean and the median of the distribution. In 
practice, for our final values of the effective decay rates (shaded entry in 
Tab. |2|), we look for the "best fits", i.e. the ones closest to the mean and 
median of the weighted distribution and satisfying the constraint that their 
confidence level is at least 80%. The statistical error is then the jackknife er- 
ror of the best fit, while the systematic error is given by the 68.3% confidence 
interval of the weighted distribution. The final error is computed by adding 
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both errors in quadrature. We also analyzed the effect of adding a cut in the 
x axis in order to avoid possible finite volume effects at large values of the 
taxi driver distance (cut at ||x|| = L or :y y^, i.e. the maximum distance in 
a L 3 x (2L) box in the continuum, with periodic boundary conditions), and 
we found results that are fully compatible with the ones without the cut, 
however, with a decrease of the systematic error. 

Combining all these possibilities to extract the effective decay rate by 
using the mean and median from the weighted distribution (labeled mean 
and median in Tab. [2]) using either the confidence level or x 2 /d.o.f., mean 
and median from "best fit" labeled fit(mean) and fit(median) in Tab. [2]), we 
thus have in total 8 different ways of computing the effective decay rate and 
its statistical @ and systematic error, see Tab. [2J 

In Fig. [21 we show the comparison of the final results for the effective decay 
rate obtained using these 8 different methods, for the case of the ensemble 
Eg. All methods lead to compatible results both in the central value, as 
well as in the total error, which gives us confidence that we indeed have the 
systematic errors of the effective decay rate under control. 

For the continuum limit scaling test of the effective decay rate, we decided 
to use method 4 of Tab. [21 taking the value for p that yields a fitted p as 
close as possible to the median of the CL-weighted histogram. However, we 
want to emphasize that all other variants would be equally good choices. 

Before discussing the numerical results for the effective decay rate, we 
show in Fig. [3] the exponential decay of the maximum norm of the overlap 
Dirac operator as a function of the taxi-driver distance. We use four different 
lattice spacings, ranging from 0.079 to 0.042 fm, while keeping the value of 
L/a = 32 fixed, as discussed above. We observe that decreasing the lattice 
spacing leads to an increasing slope of ||-D v|| max - This signals the continuum 
limit restoration of point-like locality. 

In Tab. [31 we summarize the results obtained for the values of the decay 
rate p, together with their statistical, systematic and total errors. As already 
suspected from Fig. [I] the total error is dominated by the systematic part. 

Using data from Tab. [31 we perform the continuum limit scaling test of 
the inverse decay rate (locality radius) 1/p. The values at non-zero lattice 
spacings can be interpreted as the physical length that correspond to the 
slowest decrease of the Dirac operator norm. The value of the locality radius 
extrapolated to a = is 0.007(9) fm, i.e. compatible with zero. Thus, in 
the continuum limit, point-like locality is indeed restored. If we consider the 

4 The statistical error quoted for the values of the mean and the median extracted from 
the distributions were obtained by analyzing the weighted distribution of all the statistical 
errors of the different fits computed with jackknifc. 
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Table 2: Different methods used to 
calculate the value of the effective 
decay rate p and the corresponding 
error. 



> 

CD 



le-05 



== le-10 



le-15 



le-20 




Figure 3: Exponential decay of the maximum norm of the overlap Dirac 
operator for ensembles -8^32, Ce,32, -D/1,32, Eg. 

lattice artefacts and therefore perform a fit adding a term quadratic in a, 
instead of a simple linear fit, we obtain a compatible result: -0.002(37) fm, 
but with a much larger error, indicating that we are not sensitive to 0(a 2 ) 
cut-off effects in the effective decay rates. 
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[2.923,3.634] 


W,32 
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0.333 
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[1.953,2.205] 


Dh,32 
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0.094(1) 


[1.326,1.785] 


Ei 


12.677 


0.056 


0.311 


0.316 


0.079(2) 


[1.176,1.302] 



Table 3: Results for the decay rate p and the corresponding errors: statistical 
<5 stat p, systematic S syst p, total 5p. We also give the inverse decay rate (with its 
combined statistical and systematic error in parentheses) and the fit interval 
for the best fit. 
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Figure 4: Continuum limit scaling of the inverse decay rate 1/p. 

3.3 Comparison of decay rates and hadron masses 

The relevance of the decay rate p of the overlap operator can be quantified 
when compared to the size of a hadron. The locality radius, r\ oc = 1/p, is 
expected to be sufficiently smaller than the hadron size in order to extract the 
hadron mass without being affected by non-locality problems. To illustrate 
this point, we consider the approximation in which the size of the hadron is 
related to the inverse of its mass m. There is a similarity of the asymptotic 
exponential decay of zero-momentum correlation functions (with Euclidean 
time) and the likewise asymptotic exponential decay of the overlap Dirac 
operator norm (with taxi-driver distance). At large Euclidean time or taxi- 
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driver distance one finds (barring effects from periodic boundary conditions): 




C(t,p = 0) = Ae~ 




(17) 



where m corresponds to the mass of the lightest hadron with given quantum 
numbers and p e g denotes the effective decay rate. 

In this simplified picture, in order to have an extraction of a hadron 
mass which is not influenced by a possibly too slow decay rate of the overlap 
operator, the decay rate is expected to be larger than the measured hadron 
mass: p > ma. 

In terms of the Compton wavelength of the hadron Ac = 1/m and the 
locality radius r\ oc = 1/pwe want Ac > Hoc- When the condition Ac > ri oc is 
met, once we have reached the asymptotic Euclidean time region to extract 
masses, the norm of the overlap operator has decayed already sufficiently 
strongly such as not to influence the hadron mass measurement. 

We now show that the condition Ac > r\ oc is indeed satisfied in our 
mixed action simulations, considering the masses of the lightest particles in 
the meson and baryon sectors. We take the unitary values of the pion and 
nucleon mass, given in Ref. [53]. Since we match the pion mass in our mixed 
action setup, the values of the pion mass are the same as in the unitary setup, 
by definition, whereas the masses of the nucleon can differ from the unitary 
ones by 0(a 2 ) effects. These turn out to be very small in practice - see Sec. O 
for details. We summarize the values for the pion mass, the nucleon mass, 
the decay rate and the ratios of the masses and the decay rate in Tab. 21 

In Fig. [51 we show the ratio am/ p as a function of the lattice spacing 
a in a fixed physical situation, see above. Both for the cases of the pion 
and the nucleon masses, the values of the ratio are below one (slightly above 
one in the case of the nucleon at = 3.9) and the continuum limit value 
is compatible with zero. This suggests that the overlap Dirac operator is 
sufficiently local in our mixed action setup to allow for a clean computation 
of these masses. However, working at a coarser lattice spacing or considering 
heavier hadrons, one could enter the regime where the size of the hadron of 
interest is smaller than the locality radius. When using overlap fermions, 
it is thus crucial to monitor that the locality radius is sufficiently smaller 
than the hadron size for all the values of the lattice spacing which are being 
considered. 
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am N /p 
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0.159(2) 


0.283(24) 


0.511(6) 


1.056(90) 


Ct 


0.121(4) 


0.221(11) 


0.409(6) 


0.908(41) 


D e 


0.098(2) 


0.182(4) 


0.305(4) 


0.756(28) 


E( 


0.075(2) 


0.141(5) 






cont. limit 




-0.025(21) 




0.17(16) 


cont. limit (quadratic) 




-0.10(9) 




-0.3(1.2) 



Table 4: The pion masses m n and nucleon masses in lattice units and 
the ratios of these masses divided by the decay rate of the overlap Dirac 
operator am n /p and amjv / 'p. The continuum limit values quoted in the table 
correspond to a linear extrapolation in the first case and in the second case 
we added a quadratic term to take into account possible higher order lattice 
artefacts. 
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Figure 5: Continuum limit scaling of am/ p for the pion (left) and nucleon 
(right) masses. In these plots only the linear fits are shown. 

3.4 Dependence on the s parameter 

The overlap Dirac operator is constructed from the Wilson Dirac operator 
taken at large negative mass shift, depending on the s parameter (see Eq. (j3J), 
which has to fulfill \s\ < 1, as shown in Ref. |54] . In this subsection, we 
analyze the influence of the s-parameter on the decay rate p of the overlap 
operator norm. This allows to choose the optimal value of s for simulations, 
i.e. the one that gives the highest value of p. We compare the case when 
our gauge field configurations were HYP smeared (1 iteration) - which is the 
case for the results in the previous sections - with the one of un-smeared 
configurations. Since HYP smearing brings the plaquette values closer to 
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Figure 6: Exponential decay of the norm of the overlap Dirac operator 
(normalized by the value at ||x|| = 0) for several values of the s parameter 
as a function of the taxi-driver distance. The graph represents our results 
for the Ensemble Be. (left) 1 iteration of HYP smearing was applied to our 
gauge field configurations, (right) No HYP smearing applied. 

unity and lifts the low-lying eigenvalues of the kernel Wilson-Dirac operator, 
improved locality can be expected and is one of the motivations of using HYP 
smearing in our mixed action setup. 

In Fig. [HI we plot the dependence of the norm of the overlap Dirac operator 
on the taxi-driver distance, for the ensemble with the coarsest lattice spacing 
Be. We compare several values of the s parameter for HYP smeared and 
original, non-smeared configurations. In all cases, we observe an exponential 
decay, with a significant dependence on the s parameter. The resulting values 
of the decay rate p, which are extracted using the strategy described above, 
are shown in Fig. [7J For the case of HYP smeared configurations, the optimal 
value of s is the free-field optimal value s = 0, in accordance with Ref. [55J, 
while for the original configurations s = 0.4. In general, we see a strong 
dependence of p on the s-parameter and it seems from Fig. [TJthat the main 
difference between HYP-smeared and non-smeared gauge field configurations 
is solely a shift in the s-parameter. In particular, within the errors, there 
seems to be no clear gain of using smeared or non-smeared configurations 
for the effective decay rate. Nevertheless, we do observe that in the HYP- 
smeared case the norm of the Dirac operator reaches lower values for 
the same values of the taxi driver distance than in the Non-HYP smeared 
case, as is shown in Fig. EJ 
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Figure 7: The dependence of the decay rate p on the value of the s parameter 
for HYP and non HYP smeared configurations. The lines are just to guide 
the eye. 



3.5 Finite Volume Effects 

Since our locality analysis was performed for different physical volumes (from 
ca. 1.3 to 2.5 fm), it is natural to ask whether some of these volumes are 
not too small for a reliable calculation of the decay rate. In order to see 
whether finite size effects affect the results for the decay rate, we extended 
our analysis for ensembles Bg j32 , C432, -0^,32 and Eg to the corresponding 
ensembles with smaller physical volumes: Bg 8 , Bg, Bg j2 o, Cg and D^. 

We illustrate in Fig. [S]the outcome of this analysis by comparing results 
for ensembles Cg and Cg$% (although they are practically the same at other 
values of (3). The upper plots show that in both cases we observe an ex- 
ponential decay of the overlap Dirac operator norm. However, as the lower 
plots show, for the case of L/a = 20, we observe no plateau of the effective 
decay rate and hence we cannot perform meaningful fits and extract a reliable 
value of the decay rate^|. Further investigation of the problem in free-field 
theory for several values of L/a indicates the presence of large hypercubic 
artefacts, reaching taxi-driver distances as large as ||rr|| /a ~ 20. Moreover, 
it is natural to expect that at \\x\\ / a & L/ a finite volume effects start to be- 



5 If we insist on extracting a value, it is compatible with the one of L/a = 32, but within 
a large systematic error - it therefore shows that our way of extracting the systematic error 
is reliable. We emphasize that even in the case of L/a = 32, we have to take into account 
systematic effects, which dominate over statistical errors. 
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Figure 8: (upper) Exponential decay of the norm of the overlap Dirac op- 
erator (normalized by the value at ||x|| = 0) and the effective decay rate p 
(lower) as a function of the taxi driver distance [fm] for f3 = 4.05, L/a = 20 
(left) and L/a = 32 (right). In the lower left plot, we observe no plateau. In 
the right plot, the plateau is present and hence we can perform meaningful 
fits. 



come large. Hence, the plateau region of the effective decay rates is located 
approximately between ||x|| /a = 20 and 32, as observed in the lower right 
plot of Fig. [HJ This conclusion is valid for all values of (3, as can be seen 
in Fig. [T]and it motivated our choice of L/a = 32 for all lattice spacings, 
rather than keeping the physical extent L fixed. Indeed, we observe that 
the decay rate does not depend in a significant way on the volume once the 
lattice size is sufficiently large. Therefore, for the range of lattice spacings 
considered in this work, we expect that the finite size effects in the decay 
rate are negligible when working on L/a = 32 lattices. We have also checked 
that performing the continuum limit with the physical volume fixed to 1.3 
fm yields a compatible result, but with much worse precision, due to large 
hypercubic artefacts for L/a = 16 (/3 = 3.9) and L/a = 20 (P = 4.05). 
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4 Pion decay constant and the role of the zero 
modes 



In Ref. [27J, we have reported the results of a continuum limit scaling test 
of the pion decay constant f n , using maximally twisted mass (MTM) sea 
fermions and two kinds of valence quarks: MTM (unitary setup) and overlap 
(mixed action setup). In the mixed action setup, both kinds of fermion ac- 
tions were matched, using the pion mass or by employing alternative match- 
ing conditions. While in the matched mixed action setup the value of 
(or any other observable) at finite lattice spacing might be different for two 
different fermion discretizations, one expects that this difference vanishes 
in the continuum limit, provided that the matching conditions were chosen 
properly. Let us shortly summarize the findings of Ref. |27| . 

In order to compare the continuum limit of f n for the unitary and the 
mixed action setups, we matched them by employing several matching con- 
ditions. 

1. Naive matching condition - we extracted the pion mass from the cor- 
relator constructed from the pseudoscalar interpolating field (PP cor- 
relator Cpp(t)) for both MTM and overlap valence quarks. We defined 
the matching point as the overlap valence quark mass that yields the 
same PP correlator pion mass as in the unitary case. 

2. Improved matching condition - the unitary pion mass was extracted 
from the PP correlator, but the overlap valence pion mass was ex- 
tracted from a correlator defined as Cpp-ss{t) = Cpp(t) — Css(t). In 
such a correlator (PP-SS correlator), the effects of zero modes exactly 
cancel, since zero modes couple in the same way to the pseudoscalar 
and scalar correlators. We defined the matching point as the overlap 
valence quark mass that yields the same PP-SS correlator pion mass 
as the PP correlator pion mass of the unitary case. 

3. Alternative matching condition - matching was performed using the 
PP correlator, but at a heavier quark mass, using partially quenched 
MTM data. In this way, one also suppresses the effects of zero modes, 
since their leading contribution is proportional to the squared inverse 
quark mass. 

For the regime of parameters that we had chosen, we found that using the 
naive matching condition leads to incompatible continuum limits for r /* m 
and ro/° v , if both are extracted from the PP correlator. This was attributed 
to the effects of zero modes that couple strongly to this correlator when 
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Figure 9: Continuum limit scaling of rof n . We used the improved match- 
ing condition and extracted /° v from the PP-SS correlator. The solid line 
represents a linear extrapolation in a 2 to the continuum limit. 
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0.0749(16) 
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0.004 





Table 5: The pseudoscalar masses and overlap matching quark masses for 
different ensembles and different matching conditions. Results for ensembles 
Bg, Cg and Dg are taken from Ref. [27J. The ensemble Eg at a very fine value 
of the lattice spacing has been newly analyzed in this work. 

overlap valence quarks are used. These zero mode contributions are not 
matched by the MTM sea quark action, at least not at the values of the 
lattice spacing used in Ref. [27]. To test this hypothesis, we employed the 
improved matching condition and extracted /° v from the PP-SS correlator. 
In this way, we found compatible continuum limits. By investigating the 
three matching conditions described above, we finally concluded that each of 
them can lead to compatible continuum limits of /* m and f° v , provided that 
the latter is extracted from the PP-SS correlator, which does not have the 
contribution of zero modes. 

The aim of the remainder of this subsection is to update this analysis 
with the fourth lattice spacing (ensemble Eg at j3 — 4.35). We followed the 
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same procedure as in the case of the coarser lattice spacing ensembles (Bi, 
Ce, Di) - we extracted the quark mass dependence of the overlap pion mass 
from the PP-SS correlator and found the matching mass, employing thus the 
improved matching condition. The values of the matching mass, employing 
different matching conditions are summarized in Tab. [5j At this matching 
mass, we computed the pion decay constant f° v from Eq. and compared 
it to the unitary value /* m , computed from Eq. ([S]). The results for our 4 
lattice spacings are presented in Fig. [9j 

We obtained a result which is fully consistent with our finding in Ref . [27] . 
Both the unitary setup and the mixed action setup lead to compatible con- 
tinuum limits (compatible between themselves and with the values found in 
Ref. [27 \). Let us emphasize again that the continuum limits are compatible 
only when /° v is extracted from the PP-SS correlator, i.e. when the contri- 
bution of zero modes is removed. The necessity of using proper matching 
conditions as found in our work should, in our opinion, be taken into ac- 
count for any calculation in a mixed action setup which uses chiral invariant 
fermions in the quark sector. 

In order to further corroborate this conclusion we also tried the naive 
matching condition, using solely the PP correlator, for a larger volume with 
L pa 2 fm and a pion mass of w 480 MeV, i.e. the ensemble B s of 
Tab. [TJ For this ensemble, we could observe a significant drop of the difference 
between f° v and /* m from about 40% for the small volume and small pion 
mass to about only 8% for the ensemble B s . Details of this analysis are 
presented in Appendix IA1 

5 Low energy constants of mixed action chiral 
perturbation theory 

5.1 Pseudoscalar meson masses and A Mix 

The chiral Lagrangian including 0(a 2 ) corrections for a mixed action with 
twisted mass sea and Ginsparg- Wilson valence quarks can be derived from 
the twisted mass chiral Lagrangian |56[ [57J [581 ESI [60] . In this section we will 
use the expressions of Ref. [61]. They contain the Wilson chiral perturbation 
theory (%PT) low energy constants (LECs) Wo, Wq 7 8 as well as an additional 
LEC, Wm-i appearing in a mixed action. A dependence on the twist angle 
uJq, defined by atan(wo) = fi/m is also present. 

Considering the power-counting a 2 ~ m leads to the following leading 
order (LO) expressions for sea and valence quark mass dependence of the 
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light pseudoscalar meson masses (521 E31 EH ED E51 EH] (we set u = 7r/2): 



Mi s , ± = 2B 0f i, (18) 

32 

^s 2 s,o,co„ n = 2B^-a 2 -W^ (19) 

M vv = 2B m OY , (20) 

M vs = B (m ov + fi) + a 2 ^W M -a 2 ^W^ (21) 

where the convention for the pion decay constant in the chiral limit / is 
such that f n is 132 MeV, B is a low energy constant related to the chiral 
condensate: B = —2{0\uu\0)/f 2 |67j, a = 2W a, Mss,± is the charged sea- 
sea (SS) meson mass, Mss,o,conn the neutral sea-sea meson mass computed 
solely from the quark-connected contributions and Mvv, Mvs are valence- 
valence (VV) and valence-sea (VS) meson masses, respectively. 

Rearranging the above expressions, one can find the relation between the 
dimensionless combination of LECs {tqWq)(Wm ~ 2W%) and the SS, VV and 
VS meson masses: 

{r°Wo ) {W M ~ 2W 8 ) = r I M vs - J [~ ) ■ (22) 

This equation will allow us to determine the combination {tqWq)(Wm — 2W£) 
from our lattice data for the masses M V v, Mss,± an d M V s- 

An alternative parametrization of the VS meson mass [621 |63j El] is: 

M vs = B (m ov + a 2 A Mix . (23) 

The relation between (Wm — 2WQ and Amix is the following: 

A-Mix — ~p ' V ' 

For our ensembles B^, Cg, and Bh, we have computed the pseudoscalar 
correlation functions for mesons constructed from two MTM quarks (SS,±), 
two overlap quarks (VV) and a combination of one MTM and one overlap 
quark (VS). We will call the latter mixed correlators. From each of these 
correlators, we have extracted pseudoscalar meson masses. The dependence 
of these masses (Mgs,±, Mvv and Mvs, respectively) on the quark mass 
is shown in Fig. [TOj Since the combination of LECs (^W^Wm — 2Wg) 
and hence also Auix are determined by the difference between M vs and the 
average of M vv and M| s± , we use the so determined meson masses in the 
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Figure 10: The quark mass dependence of the mixed (VS) pseudoscalar me- 
son mass and the overlap (VV) pseudoscalar meson mass for the ensembles of 
Tab. [TJ The latter has been shifted by the unitary MTM (SS,±) pseudoscalar 
meson mass, which does not depend on m ov . The quark mass independent 
vertical distance between 2(aMvs) 2 and (aMw) 2 + (<^fss,±) 2 determines the 
value of A Mix, which should in the LO of considered here not depend 
on the quark mass. 



graph. At the considered order of xPT, this difference should not depend on 
the overlap quark mass m ov and the twisted quark mass /i, i.e. the data for 
My S and My V + M| s ± should be parallel as a function of the quark mass. 

As the Fig. [10] shows, this is indeed the case for all ensembles we have 
investigated. In Tab. [6J we also show the extracted values of the combination 
( r o^o )(Wm— 2Wg) and the corresponding values of Amix- The overlap quark 
mass dependence of these LECs is compatible with a constant behaviour and 
hence we have used constant fits of the difference 2My S — (M vv + M| s± ) 
to extract {tqWq){Wm — 2Wg) for each ensemble, see also Fig. [TTJ We find 
that the extracted values are compatible for all ensembles, as they should at 
his order of x?T. In fact, the spread in {tqWq)(Wm — 2Wg) (or in Amix) 
between our light sea quark ensembles, labeled with the subscript £, is very 
small, suggesting that lattice spacing effects are rather small, at least within 
the precision we could reach here. 

As our final values for (tqW^CWm— 2Wg) and A M ^ X , we take the weighted 
averages of the values from Tab. |6] (with the inverse statistical error squared 
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Figure 11: The quark mass dependence of the combination of LECs 
{tqWq){Wm ~ 2H^g) for the ensembles of Tab. [TJ The vertical scale is the 
same on all plots, indicating that the values extracted for all ensembles are 
compatible with each other. 



Ensemble 


matching 

arriov 


(r 6 <) (W M 
at matching mass 


-2Wi) 
constant fit 


<-^Mix 

[MeV] 


Be 


0.007 


0.092(23) 


0.081(12) 


997(37)(21) 


Ci 


0.005 


0.075(22) 


0.068(11) 


951(38)(21) 


D £ 


0.002 


0.081(33) 


0.079(13) 


968(40)(21) 


B h 


0.015 


0.056(11) 


0.054(9) 


901(38)(20) 



Table 6: The extracted values of the combination of LECs 
{tq Wq ) (Wm — 2Wg) . We show the values at the matching point 
Myv = ^fss,± (the matching quark mass is given in the column am ov ) and 
from the constant fit to 2(aM V s) 2 — (aM V v) 2 + (aMss,±) 2 for the whole 
quark mass range shown in Fig. [TTJ The last column is Aj^j, in MeV, for 
which the first error is statistical and the second one the systematic error 
of conversion to MeV, coming from uncertainties in the values of lattice 
spacings. 
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as the weight). This yields: 



(r 6 W*){W M - 2W' & ) = 0.067(6)(14), 

A Afe = [951(21)(50)MeV] 4 , 

where the first error is statistical (error of the weighted average) and the 
second one is systematic, which we take as one half of the spread between 
the values for the ensembles Bi and B^. This implies an approx. 140 MeV 
splitting between the VS and VV meson masses at the matching mass, for 
ensemble B^. Note that our result for Wm ~ 2Wg is compatible with the 
constraint W M - > 0, found in Ref. [65] . 

The value of Amix that we have determined is comparable to, but some- 
what larger in comparison with previous studies with different mixed action 
setups: A Mix = [706(4) MeV] 4 [68] and [678(13) MeV] 4 (with a finer lattice 
spacing) |21| for domain wall fermions on a staggered sea, [769(77) MeV] 4 and 
[861(90) MeV] 4 (for a heavier pion mass) [T7] for overlap fermions on a clover 
sea and -[427(338) MeV] 4 [Tg] or [416(27) MeV] 4 [69] for overlap fermions on 
a domain wall seal. 



5.2 Unitarity violations in the scalar correlator 
5.2.1 Theoretical predictions 

Mixed actions, even when related to each other by suitable matching con- 
ditions, violate unitarity at any non-vanishing value of the lattice spacing. 
Unitarity is then only restored in the continuum limit. So far, we have looked 
at unitarity violations in the pion decay constant at the matching mass and 
in mixed correlators. Here, we investigate yet another way of studying uni- 
tarity violations, through a quantity which is known to depend strongly on 
them and, in principle, allows to quantify their size - the flavour non-singlet 
scalar correlator. The reason is that within x?T, unitarity violations appear 
as unphysical double poles in the propagator, affecting strongly the flavour 
non-singlet scalar correlator even at the matching point, i.e. if the valence- 
valence (VV) pseudoscalar meson mass is matched to the sea-sea (SS) meson 
mass [TOj EJ UM E3]. In fact, the scalar correlator can turn negative in a 
mixed action setup which is clearly an unphysical behaviour. 

6 In Ref. [SS], the constraint takes the form 2Wm — Wg > 0. However, Wm in our 
notation receives a factor of 4 compared to the one defined in Ref. [65J . 

7 The former value has been obtained using a different (indirect) method - by examining 
a meson state that wraps around the time boundary, instead of a direct extraction from 
mixed correlators. 
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The flavour non-singlet scalar correlator can be approximated by: 

C™(t) = -2 (V S 2 S;± - M vv + a 2 ^ W^j ■ B DP (T, L, t, M vv , M vv ) 

-2B SP (T, L, t, M vv , M vv ) + 2B SP (T, L, t, M vs , M vs ) 
+Ae -ma °* + excited states , (25) 

where Ssp and Sdp are the single and double pole bubble functions including 
finite volume effects |6T] and the last terms describe contributions from scalar 
mesons (the lightest of which is ao ). The label VV in C™(t) reminds that 
overlap valence quarks are used in the scalar correlator. 

Assuming matching of pseudoscalar meson masses Myy = Mss,±, in the 
limit T — > oo, and at large Euclidean times t, the bubble contributions can 
be simplified to yield: 



B 2 

^sca (t) ~~ ^ -^jj 



e 



-2M vs t e -2M vv t / ]_g 



M vs M vv 



16 \ 

M 2 y + a 2 - W^l+Myyt)) 



.(26) 



Identifying the LEC 7ss in Ref. [73] with 7ss = y| WoWg, one recovers, in 
the large t limit, the corresponding expression of Eq. (23) of Ref. [75] . 



5.2.2 Lattice results 

Due to the negative coefficients appearing in (j26|) it can be that the SS 
correlator becomes negative at large Euclidean times thus indicating a lack 
of unitarity. However, in order to isolate the different contributions to the 
violations of unitarity, it is useful to disentangle the effect of zero modes, 
which can have a potentially very large contribution to the SS correlator. 
This contribution is actually equal to the zero modes contribution to the 
PP correlator. As previously discussed, this is behind the choice of the 
improved matching condition, which allowed us to extract the pseudoscalar 
meson mass and decay constant by computing the PP-SS correlator, where 
the effect of zero modes is exactly canceled. In Fig. [121 we show the scalar 
correlator for ensembles Be and Df . We have computed the index of the 
overlap Dirac operator for each configuration and we have calculated the SS 
correlator for all configurations with a given index. For large Euclidean times 
(i.e. in the region where Eq. (j26|) should be applicable), the SS correlator 
has a tendency to become more and more negative for configurations with 
an increasing number of zero modes, i.e. in higher topological charge sectors. 
When the correlator is evaluated on all configurations, the SS correlator is 
clearly negative, which signals unitarity violations. 
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Figure 12: The scalar correlator C^(t) in fixed topological charge sectors 
(index 0,1,2,3,4) and from summing over all topological charge sectors (all). 
In addition, we show the results for the scalar correlator with explicitly sub- 
tracted zero modes (SS sub.). Upper: ensemble Bi, (3 = 3.9, L/a = 16, 
am s = 0.004, am ov = 0.009. Lower: ensemble D h (3 = 4.2, L/a = 24, 
am s = 0.002, am ov = 0.004. 
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Figure 13: Fits of Eq. (j25p . using ensemble B s , for which the zero modes 
effects are very much suppressed. Blue solid line is the full fit of Eq. (|25|) . 
which is a sum of 4 contributions: double pole (DP) - black solid line, single 
pole containing Myy (SP VV) - gray dashed line, single pole containing My$ 
(SP VS) - cyan dash-dotted line, scalar a meson (aO) - red dotted line. 

In order to determine the low energy constant W% appearing in Eq. (|25|) . 
it is much safer to use larger volumes and quark masses, i.e. a situation when 
the contribution of zero modes is highly suppressed. Among our ensembles, 
B s has the desired properties and will be used to extract W%. 

Another possibility is to restrict oneself to topologically trivial configura- 
tions. However, it is clear that the number of such configurations is rather 
small for each ensemble (between 10 and 20% of all configurations) and hence 
the signal is of poor quality, as can also be inferred from Fig. [T2] directly. 
Moreover, taking only topologically trivial configurations means that the re- 
sults of the ensemble average will differ from the full result (i.e. from all 
topological sectors) by power-like finite volume effects [35l ITS] . Therefore, 
such an analysis will be used only as a cross-check and as an estimate of a 
systematic error related to our computation using ensemble B s . 

Yet another way to proceed is to explicitly subtract the contribution 
of zero modes at the level of propagators. This is a potentially dangerous 
procedure with hard-to-control systematic effects. The SS correlator with 
explicitly subtracted zero modes (which we call SS subtr.) is shown in Fig. [T2J 
However, given the fact that the zero mode subtraction is a rather doubtful 
procedure, we will quote the value obtained with this method, but we will 
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not take it into account in our analysis leading to the final value of W$. 

We now describe our strategy to extract Wg from ensemble B s , using 
ensembles Bg, Bh, Ce and Dg to estimate the systematic error. We fit Eq. (1251) . 
which has 3 fitting parameters: the LEC of interest - Wg, the amplitude of 
ao scalar meson contribution and its mass m ao (however, we do not attempt 
to extract the ao contribution quantitatively). In addition, we add a fourth 
fitting parameter - the overall normalization of the correlator N, into which 
we absorb unknown renormalization factors Z$- We work with symmetrized 
correlators and take into account the finite extent of the system by adding 
in Eq. (|25() the terms with t — > T — t. 

The results of our fits for ensemble B s are shown in Fig. [I3J The di- 
mensionless combination of LECs rgW^VKg determined from our best fit is 
—0.00642(15), where the error is statistical only. We now want to deter- 
mine the systematic error, related to: residual discretization effects, finite 
volume effects, choice of the Euclidean time fit range and finite precision of 
the input parameters of our fits: aMyy = 0.1918(12) at the matching mass 
am ov = 0.017 (see Appendix Ej), aM ss ,± = 0.19403(50)0 [B], A Mix (taken 

8 Note that aMyy is not exactly equal to aAfgg,±, which we take into account by 
including the difference M| s ± — My V in the fits of Eq. (f25j) . 
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central value 


-0.0064 


B s 


fit range 


0.0015 


B s 




0.0004 


B s 


aM vv 


0.0002 


B s 


aM ss 


0.0001 


B s 


including/excluding Mf s ± — My V 


0.0003 


B s 


discretization effects 


0.0008 


Be, Ce, De 


finite volume 


0.0009 


B s , B h 


sea quark mass 


0.0012 


Be, B h 


ro/a, r f 


0.0006 


Ref. p] 


total systematic 


0.0024 


B s , B h , Be, Ce, De 


statistical 


0.0002 


B s 



Table 7: Error budget for our estimates of the combination of LECs rgW^ Wj^. 
The last column shows which ensembles were used to determine the given 
type of systematic error (or central value). The total systematic error comes 
from combining the individual ones in quadrature. 

from the previous subsection), r /a = 5.25(2), r f = 0.259(11) |39| . 

We address the errors related to discretization effects, finite volume effects 
and sea quark mass using our small-volume ensembles at 3 lattice spacings: 
Be, Ce, De (at a pion mass of about 300 MeV) and Bh (approximately 450 
MeV). We illustrate fits for these ensembles in Fig.[TH The remaining system- 
atic errors were estimated using ensemble B s only and varying the respective 
quantities, in particular the fitting range. We have also checked the stability 
of our results with respect to including or excluding the 2 fitting parame- 
ters related to the scalar meson (ao, A). This gives fully compatible results 
(within the systematic error from the choice of the fitting range), however 
we are then limited to larger values of t/a. The summary of all systematic 
errors is presented in Tab. [7J 

The final value that we quote is: 

rX^s = -0.0064(2) (24), 

where the first error is statistical and the second systematic. This value can 
be compared to the recent determination (with twisted mass fermions) [75] 
from the connected contribution to the neutral pion correlator (Eq. (|T9l ). 
which gives ^o^o^s = — 0.0106(11) (ensemble Be at a larger volume with 
L/a = 24 instead of 16) or r^W 2 Wg = -0.0096(22) (ensemble B e ). We also 
give here as a cross-check the value obtained from explicit subtraction of 



32 



(jW!) (W m - ZW[) 0.067(6)(14) I A Mu , [951(21)(50) MeV] 4 
(r$Wjf)W£ -0.0064(2) (24) w' s -[529(4) (51) MeV] 4 

(r%W$)W M 0.054(6)(15) w M [901 (21) (62) MeV] 4 



Table 8: The values for different combinations of LECs determined in this 
work. The first error is statistical and the second one systematic. 

zero modes: ^o^o^s = — 0.0127(8) (where the error is purely statistical). 
Although this value has the right order of magnitude, we emphasize again 
that the validity of the explicit zero modes subtraction procedure is doubt- 
ful and hence this result should be interpreted with caution. We also note 
that our finding for W$ is compatible with the constraint W$ < found in 
Refs. [661 [761 [771 [78]. 

Our results for the combination of LECs ToW^ 2 W 7 ^ and tIWq(W m — 2Wg) 
can be combined to isolate: 

r 6 W*W M = 0.054(6)(15), 

where the meaning of the two errors is as above. 
We also define quantities analogous to A Mix : 

< - ™, (27) 

«* = i!». (28) 

With such definitions, the following relation holds: Amix = Wm — 2«7 8 . Our 
results imply the values: w M = [901 (21) (62) MeV] 4 and w' 8 = -[529(4)(51) MeV] 4 . 
The latter corresponds to approximately 125 MeV splitting between the 
charged and neutral (connected^) pion masses, for ensemble Bg. For conve- 
nience, we summarize all our results for mixed action x?T LECs in Tab. [HJ 



6 Light baryon masses 

In order to investigate the effects of the mixed action setup in other ob- 
servables, we have computed the nucleon and A baryon masses in the mixed 

9 If one considers the full (connected and disconnected contributions) neutral pion mass 
Mss,o, the splitting Afss,o — Mss.± is governed by a combination of LECs W$ and W' % (see 
e.g. Ref. [60J). The LEC Wq is not directly accessible from the overlap scalar correlator. 
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action and the unitary setup. Given the numerical cost of the inversion of the 
overlap operator, we performed the analysis on the following subset of gauge 
ensembles generated on only small lattice sizes up to 24 3 -48, Bp (/3 = 3.9), 
Ci {(3 = 4.05) and Dt (/3 = 4.20) that correspond to a fixed spatial extent 
of « 1.3 fm and a fixed pseudo scalar mass of ~ 300 MeV. Having three 
values of the lattice spacing available, allows us to study the continuum limit 
for the nucleon (N) and A masses in both setups. The strategy followed to 
extract the baryon masses in this section follows closely the work on the light 
|79| and strange baryon spectrum [53J by ETMC For consistency, we recall 
here the basic ingredients of the computation, the reader interested in more 
details is referred to the two aforementioned references. 

The masses of the nucleon and of A baryons are computed through the 
two-point correlators, using the following interpolating fields: 



where C is the charge conjugation matrix. The quark fields u and d refer 
to the two degenerate flavours of quarks considered in this work 1^*1 While 
the operator J N carries spin 1/2 and couples only to states which have the 
quantum numbers of the nucleon, the operator couples both to spin 1/2 
and spin 3/2 states. Note, however, that in practical lattice computations 
the spin 3/2 dominates the correlator at large time [79 J thus allowing for a 
clean extraction of the A baryon. As in nature, the spin 1/2 partner is much 
heavier than the spin 3/2 state. 

In order to improve the overlap between the ground state and the inter- 
polating fields, we employ source and sink smearing of the quark fields. We 
use Gaussian smearing |80[ I8~T] of the quarks field and APE smearing [82] 
of the gauge links entering in the Gaussian smearing procedure. We use the 
smearing parameters obtained in our previous works on baryon spectroscopy 
[T9~l 153] . Note that we use the same smearing parameters both in the twisted 
mass and overlap fermion cases. 

In order to investigate the consequences of the various matching condi- 
tions discussed in Sec. HJ we use several valence overlap quark masses. The 
baryon masses are extracted from the asymptotic Euclidean time behaviour 
of the correlators: 

10 In the unitary case, the u and d quarks refer to the fields in the so-called physical 
basis. 




e abc ( u a,T Cl ^ u c^ 
e abc ( u a,T ClfiV b) u c ? 



(29) 
(30) 
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C N > A (t) = J2^^(J NA (x,t)J N ' A (0)) - ?-^l°(J N > A (2,t-T)J N > A (0)), 

X 

(31) 

The statistical errors are estimated using 1000 bootstrap samples. Concern- 
ing the topological finite size effects, isolating the zero mode contribution 
after the fermionic integration, the leading behaviour of Eq. (l3~Tj) is: 

C(t) ~ (am ov y 3 (32) 

as discussed in Ref. [83J. 

We use the effective masses defined by: 

mf ff A (t) = - log ^^Tj = amNA + ° {e ' 5t) (33) 

where 5 is the mass difference between the ground state and the first excited 
state. In Fig. [TJ)] we show examples of effective masses of the nucleon and A 
in the mixed action setup. The results are obtained at the coarsest and the 
finest lattice spacing used in this work, corresponding to the gauge ensembles 
Be (left) and Di (right). The valence quark masses are fixed to the (improved) 
matching overlap quark masses summarized in Tab. [5] We fit the effective 
masses to a constant in the Euclidean time region where m^ A (t) become 
time independent. The corresponding fits are shown in Fig. [TS] where the 
fit result is represented by horizontal black lines and their statistical errors 
by black dotted lines. The effective mass plots in Fig. [15] show that we can 
identify a plateau region, allowing us a good determination of the nucleon 
and the A masses. 

Results for the nucleon and A masses in the unitary setup and in the 
mixed action setup (for am ov set to the matching mass) are reported in 
Tab. [9j We also show our results obtained for several valence quark masses 
in Fig. [16] for ensemble Bi (left) and for De (right) - each plot shows the 
overlap quark mass dependence of the nucleon mass (black dots) and the A 
mass (blue triangles). The matching mass is indicated by a vertical dotted 
line, while the results obtained in the unitary setup are indicated by a yellow 
band with a width corresponding to the statistical errors. Note that for very 
low valence quark masses, below the matching quark mass, a plateau region 
in the effective mass can hardly be found. This explains why statistical 
errors become larger for very small overlap quark masses. The uncertainty 
to extract reliably effective masses in this region is also responsible for the 
somewhat irregular behaviour in Fig. [16] As can be seen from the very good 
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Figure 15: Nucleoli and A effective masses in the mixed action setup. The 
valence quark mass is fixed to the (improved) matching overlap quark mass, 
see Tab. [5] for their values. In the left plot we show the results obtained for 
the ensemble (/9 = 3.9) and in the right for Di (/3 = 4.2). 



Ensemble 


N™ 




am™ 




am™ 




Be 


636 


426 


0.594(13) 


0.617(12) 


0.740(24) 


0.741(34) 


c, 


286 


142 


0.497(11) 


0.522(27) 


0.627(13) 


0.641(15) 




189 


371 


0.407(14) 


0.394(11) 


0.503(13) 


0.519(10) 



Table 9: Summary of our results for for the nucleon masses am^ and A 
masses am&. The superscript "TM" denotes the unitary setup, while "ov" 
denotes the mixed action setup. We also give the number of gauge field 
configurations iV used for the measurements on each ensemble. 

agreement of the baryon masses between the unitary and the mixed action 
setups, even in a regime where effects of chiral zero modes of the overlap 
operator are sizable in the pion sector, in the nucleon or A mass case, the 
effects are negligible for our definition of the (improved) matching mass. 

Finally, we show in Fig. [T7] the continuum limit scaling plot of the two 
baryon masses with three lattice spacings, in the unitary and mixed action 
case. All three ensembles have pion masses of « 300MeV and the same vol- 
ume of 1.3fm, so no interpolation in volume or quark masses is needed 
within our statistical accuracy to compare our results at different lattice 
spacings. For all ensembles, the results in the mixed action case and in the 
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Figure 16: Nucleoli and Delta masses vs. am ov . The vertical dotted line 
indicates the matching mass. The left plot shows results for ensemble Bi 
(/3 = 3.9) and the right one for D e (/3 = 4.2). 



unitary case are compatible within errors in the continuum limit. Further- 
more, in the two setups, the results for the nucleon and for the A mass 
show only a small dependence on the lattice spacings. We perform a linear 
fit in a 2 represented by dotted lines (red in the unitary case and black in 
the mixed action one). The continuum extrapolated values are represented 
slightly shifted for better readability. 
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Figure 17: The continuum limit scaling of the nucleon and A masses. 
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7 Summary and prospects 



In this paper, we have studied a particular example of a mixed action setup, 
namely valence overlap fermions on maximally twisted mass sea fermions. 
By adding a fourth and finer value of the lattice spacing than available to us 
earlier [27], we were able to perform a continuum limit scaling test for the 
pion decay constant and to compare the unitary and mixed action theories. 
We found a continuum limit value of the pion decay decay constant that is in 
full agreement between the two theories and confirmed thus the conclusion 
of Ref. (27J that with an improved matching condition the dangerous effects 
arising from the mismatch of the chiral zero modes - between the chiral 
invariant overlap Dirac operator in the valence sector and the twisted mass 
operator in the sea - can be avoided 

Applying the same strategy for other observables for which we took the 
nucleon and A masses, we found that also these quantities are not affected 
by the presence of the zero modes, once the improved matching condition 
is applied. This supports the interpretation reached in Ref. [27J that the 
improved matching condition leads to a proper matching of the sea and 
valence theories. In general, in the context of mixed actions, we nevertheless 
advocate to check the possible effects on physical observables of the mismatch 
of zero modes between valence and sea sectors, introduced by the use of a 
valence overlap Dirac operator. 

We also performed a continuum limit investigation of the locality proper- 
ties of the overlap Dirac operator. While at any non-zero value of the lattice 
spacing the overlap Dirac operator is exponentially localized, our study, em- 
ploying four values of the lattice spacing, strongly suggests that a point-like 
localization is indeed recovered in the continuum limit. One essential ele- 
ment in our investigation was to compare the exponential decay rate of the 
norm of the overlap Dirac operator with hadronic scales, for which we took 
the pion and nucleon masses. We could demonstrate that for all values of 
the lattice spacing, the ratio of the hadronic masses over the decay rate is 
smaller than one and assumes a value compatible with zero in the continuum 
limit. Thus, the exponential locality of the overlap Dirac operator is not 
expected to distort the evaluation of these hadronic observables. We think 
that such kind of tests should be performed on hadronic physical quantities 
one is interested in, in order to monitor the possible distortions of the results 
by a finite decay rate of the overlap Dirac operator. 

Finally, we compared our numerical computations of meson masses with 
predictions of mixed action %PT and extracted in this way some of the LECs 
which parametrize the mixed action chiral effective Lagrangian. We also 
provided evidence that there are indeed unitarity violations in the mixed 
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action setup and we quantified these effects by computing the Amzx and W$ 
parameters. 

In summary, we believe that the results of Ref. \27\ . in combination with 
the investigations performed in this work, can serve as a basis for future mixed 
action computations employing chiral invariant fermions in the valence sector: 
we advocate to use the improved matching condition proposed in Ref. [27J; we 
suggest to test whether the localization rate of the overlap Dirac operator is 
sufficiently larger than the hadronic scale of interest; we further point to the 
possibility to use mixed action chiral perturbation theory to look at effects of 
a mixed action setup, allowing in particular to isolate the effects of unitarity 
violations. 

Acknowledgments We thank the European Twisted Mass Collaboration 
for generating ensembles of gauge field configurations that we have used for 
this work and for a very enjoyable collaboration. We are grateful to 0. Bar 
for useful comments and suggestions. Sec. [3] of this work originated from a 
DESY summer student project carried out by 0. Haas, E. Panzer, S. Tord- 
jman. The computer time for this project was made available to us by the 
Leibniz Rechenzentrum in Munich, Poznan Supercomputing and Network- 
ing Center in Poznah, GENCI-CINES grant 2010-052271 and CCIN2P3 in 
Lyon. We thank these computer centers and their staff for all technical ad- 
vice and help. K.C. has been supported by Foundation for Polish Science 
fellowship "Kolumb" and Ministry of Science and Higher Education grant nr. 
N N202 237437. This work has been supported in part by the DFG Son- 
derforschungsbereich/Transregio SFB/TR9. G.H. acknowledges the support 
from the Spanish Ministry for Education and Science project FPA2009-09017, 
the Consolider-Ingenio 2010 Programme CPAN (CSD2007-00042), the Co- 
munidad Autonoma de Madrid (HEPHACOS P-ESP-00346 and HEPHACOS 
S2009/ESP-1473) and the European project STRONGnet (PITN-GA-2009- 
238353). 

A Results for ensemble B s 

One of the main conclusions of Ref. |27| were the values of the relevant 
parameters - physical volume and pion mass - for which the role of the 
zero modes is strongly suppressed. With such parameters, the simulation 
results are safe against the contribution of zero modes, i.e. a continuum 
limit scaling test of e.g. the pion decay constant, computed with MTM 
and overlap fermions would lead to a consistent behaviour (same continuum 
limits), even when both /* m and /° v are computed from the PP correlator. 
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Figure 18: (left) Matching the overlap and MTM pion mass, (right) The 
dependence of the overlap pion decay constant on the overlap quark mass 
am ov . The matching mass is around am ov = 0.017. Ensemble B s (f3 = 3.9, 
L/a = 24, a/i = 0.0085). 

For this work, we have decided to test this expectation explicitly, em- 
ploying an ensemble with large enough volume and pion mass. We have not 
performed a full continuum limit scaling test, since this is still very expen- 
sive computationally. Instead, we compare the values of /* m and /° v at the 
matching mass, expecting only a small difference that can be attributed to 
0(a 2 ) cut-off effects. Our ensemble, labeled B s has the following parameters: 
(3 = 3.9, L/a = 24, a/i = 0.0085, which corresponds to L pa 2 fm, pa 480 
MeV in infinite volume and m^L pa 4.7. The latter is well above m v L = 4, 
the minimal product of the pion mass and spatial extent of the lattice that 
leads to negligible effects from the zero modes \27\ . 

In Fig. [181 we show the matching of the pion mass between the unitary 
setup and the mixed action setup. We find equal pion masses (m^ m = m° v ) 
at around am ov = 0.017. At this matching quark mass, we then compare the 
value of the pion decay constant. /° v is greater than /* m by 7.7% (6.8 a away). 
This difference can be compared to the differences between /* m and f° Y for 
other ensembles at the same value of the lattice spacing: 40.0% (L pa 1.3 
fm, a/i = 0.004, m^L pa 2.5), 27.0% (L pa 1.6 fm, a/i = 0.004, m^L pa 2.8), 
19.5% (L pa 2 fm, a/i = 0.004, m^L pa 3.3), 16.5% (L pa 1.3 fm, a/i = 0.0074, 
m^L pa 3.1). Indeed, for the current ensemble of interest, the difference 
of around 7.7% is much smaller than for other ensembles at j3 = 3.9, with 
either a smaller physical volume, or a smaller pion mass (or both). While 
certainly this can not be treated as a proof that the role of the zero modes 
is negligible for this ensemble, we believe that this provides a clear hint that 
our interpretation of the role of the zero modes in our overlap/MTM mixed 
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action setup is correct. 
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